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Abstract 

We present numerical simulations of vortices that appear via primary bifurcations out of the 
unstructured circular Couette flow in the Taylor-Couette system with counter-rotating as well as 
with co-rotating cylinders. The full, time dependent Navier-Stokes equations are solved with a 
combination of a finite difference and a Galerkin method for a fixed axial periodicity length of the 
vortex patterns and for a finite system of aspect ratio 12 with rigid nonrotating ends in a setup 
with radius ratio v = 0.5. Differences in structure, dynamics, symmetry properties, bifurcation and 
stability behavior between spiral vortices with azimuthal wave numbers M = ±1 and M = Taylor 
vortices are elucidated and compared in quantitative detail. Simulations in axially periodic systems 
and in finite systems with stationary rigid ends are compared with experimental spiral data. In a 
second part of the paper we determine how the above listed properties of the M = —1, 0, 1 vortex 
structures are changed by an externally imposed axial through-flow with Reynolds numbers in the 
range —40 < Re < 40. Among others we investigate when left handed or right handed spirals or 
toroidally closed vortices are preferred. 

PACS numbers: PACS number(s): 47.20.-k, 47.32.-y, 47.54.+r, 47.10.+g 
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I. INTRODUCTION 

Spiral vortices appearingin the annular gap between the concentric rotating cylinders of 
the Taylor- Couette system |l[ are a rather interesting example for the spontaneous formation 
of a helicoidal structure out of a homogeneous state of a nonlinear forced system when the 
forcing exceeds a critical threshold. Like the competing toroidally closed Taylor vortices the 
spiral vortex structures bifurcate out of the unstructured basic state of circular Couette flow 
(CCF) that is stable at small rotation rates of the inner cylinder. The spiral pattern breaks 
the rotational symmetry of the annular gap. It oscillates in time by rotating azimuthally as 
a whole thereby propagating axially. The Taylor vortex flow (TVF), on the other hand is 
rotationally symmetric and stationary. 

The spiral pattern is effectively one dimensional (like TVF) and stationary when seen 
from a co-moving frame Q|: the spiral fields do not depend on time t, axial coordinate z, 
and azimuthal angle <p separately but only via the combined phase variable cf) = kz + Mip — 
u)(k, M)t. Here k and M are the axial and azimuthal wave numbers, respectively, and u the 
frequency. In the ip — z plane of an 'unrolled' cylindrical surface the lines of constant phase, 
cf) = 0o, are straight with slope —M/k as shown in Fig.^ An azimuthal wave number M > 
implies a left handed spiral (L-SPI) while M < refer to right handed spirals (R-SPI) with 
our convention of taking k to be positive. L-SPI and R-SPI being mirror images of each 
other under the operation z — > — z are symmetry degenerate flow states. Which of them is 
realized in a particular experimental or numerical setup depends on the initial conditions. 

With the lines of constant phase in the (p — z plane being oriented for both spiral types 
obliquely to the azimuthal 'wind' of the basic CCF both spirals are advectively rotated by 
the latter like rigid objects. Their common angular velocity is <pspi — ou(k, M)/M. This 
advection enforced rigid-body rotation of the spiral vortices is also reflected by the fact 
that the axial phase velocities w p h = u/k = (pspiMjk of an L-SPI (M > 0) and of an 
R-SPI (M < 0) are opposite to each other, see Fig. ^ By the same token the rotationally 
symmetric (M = 0) structure of toroidally closed Taylor vortices is stationary {u — 0): the 
lines of constant phases being parallel to the azimuthal CCF the latter cannot advect these 
vortices. However, an externally imposed axial through-flow can advect Taylor vortices as 
well as spiral vortices. 

The external through-flow breaks the mirror symmetry between L-SPI and R-SPI. It 
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changes their rotation and propagation dynamics as well as their structural properties and 
their bifurcation behavior in different ways. This is the topic of our investigation. 

In his review [l] Tagg remarks that systematic investigation of non-axisymmetric vortex 
states that appear via primary bifurcations out of the CCF state started remarkably late 
in the history of the Taylor-Couette problem. Krueger et al. [3j predicted in 1966 primary 
transitions to non-axisymmetric rotating-wave flow which then were observed in experiments 
by Snyder ^| who had presented experimental evidence for different types of stable helical 
flow (referred to as 'spirals') a few years earlier. In 1985, an experimental survey was 
published by Andereck et al. -5] which classified a large variety of different flow states, 
including some spiral types like linear, modulated, interpenetrating, and wavy spirals etc. 
An extensive numerical linear stability analysis was then performed for a wide range of 
radius ratios by Langford et al. jf]]. At this time, Tagg et al. j^l experimentally observed 
a transition from CCF to axially standing and azimuthally traveling waves (ribbons) and 
found numerically calculated wave speeds to be in agreement with experimental results. 
Edwards [8] studied the transition from CCF to traveling waves. More recent experiments 
were done with a system of radius ratio rj = 0.5 [9J|. Spiral solutions in a relatively narrow gap 
with axially periodic boundary conditions were obtained numerically with a pseudo-spectral 
method using co-rotating helicoidal coordinates which were adapted to the expected spiral 

0. 

Various effects of an externally imposed axial through-flow in the Taylor-Couette system 
have been explored since the early 1930 so that the list of publications cannot be discussed 
here - see, e.g, Ref. [ll] for a partial and far from complete compilation. We mention here in 
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addition a few, more recent papers on selected topics beyond those listed in Ref. 
analysis of the competition between shear and centrifugal instability mechanisms 
linear SPI and TVF fronts and pulses [1^ ; weakly nonlinear bifurcation analysis of axially 
extended spiral, ribbon, and mixed vortex states with homogeneous amplitudes 

00 

; the- 
oretical/numerical investigation of the nonlinear pattern selection in the absolutely unstable 
regime under downstream evolving intensity envelopes [u]; theoretical/numerical analysis 
of noise-sustained patterns in the convectively unstable regime (related experiments are 
listed in analysis of the changes in the noise sensitivity across the convective-absolute 
stabifi* " undary 0; _; tS of ve,o Clty fields by p artl c,e image velocimetry Q. 
In this work we explore in a detailed quantitative investigation the spatio-temporal struc- 
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tures as well as the bifurcation properties of spirals and TVF in a setup with co- and counter- 
rotating cylinders of fixed radius ratio rj = 0.5 with and without an externally imposed axial 
through-flow. Most calculations were done for axially periodic boundary conditions that im- 
pose the wave length of the vortex pattern. However also a few simulations of finite systems 
with rigid, non-rotating lids were done to compare with experiments and to study the effect 
of phase propagation suppressing boundaries. The calculations were done with a time de- 
pendent finite differences method in the r — z plane combined with a spectral decomposition 
in ip which yields by construction only the stable flows. However, by selectively suppress- 
ing destabilizing modes we also were able to trace out the unstable TVF and SPI solution 
branches. We do not include in this work results on ribbons Ml, i.e., nonlinear combinations 

n 

of L and R spirals |2j] since they were unstable for the parameters investigated here. 

In Sec. |H] we present the notation for describing the Taylor Couette system and we 
describe our numerical method. In Sec. IIHI we review the spatio-temporal properties of 
TVF and SPI and we present results on their bifurcation behavior and flow structure in the 
absence of through-flow. In particular we provide detailed comparisons of the bifurcation 
and structural properties of these primary vortex states. Also comparisons with experiments 
are presented and discussed. In Sec. IIVI we elucidate the effect of an external through-flow 
on structure, dynamics, and bifurcation properties of TVF and SPI for counter-rotating 
cylinders and stationary outer cylinder. The last section contains a summary of the main 
results. 

II. SYSTEM 

We report results obtained numerically for a Taylor-Couette system with co- and counter- 
rotating cylinders. The ratio rj — T\jri of the radii r\ and r 2 of the inner and outer cylinders, 
respectively, was fixed at the value 77 = 0.5 for which also experiments have been made 
recently |9(. 

A. Theoretical description 

We consider the fluid in the annulus between the cylinders to be isothermal and incom- 
pressible with kinematic viscosity v. The gap width d = r 2 — r\ is used as the unit of 
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length and the momentum diffusion time d 2 /v radially across the gap as the time unit so 
that velocities are reduced by vjd. To characterize the driving of the system, we use the 
Reynolds numbers 

R 1 = r l VL l d/v ;R 2 = r 2 Q 2 d/u . (2.1) 

They are just the reduced azimuthal velocities of the fluid at the inner and outer cylinder, 
respectively, where Q± and Q 2 are the respective angular velocities of the cylinders. The 
inner one is always rotating counterclockwise so that Qi and Ri are positive. We explore 
positive as well as negative values of R 2 corresponding to co- as well as counter rotation 
of the cylinders, respectively. And we elucidate the effect of an externally imposed axial 
through-flow. 

Within the above described scaling, the NSE take the form 

d t u = V 2 u - (u • V)u - Vp . (2.2) 

Here p denotes the pressure reduced by pv 2 /d 2 and p is the mass density of the fluid. Using 
cylindrical coordinates, the velocity field 

u = we r + !)e fl + ioe 2 (2.3) 

is decomposed into a radial component u, an azimuthal one v , and an axial one w. 

We have solved the resulting equations subject to no slip conditions at the cylinders. In 
Sec. IIIIDI we present simulations of systems with axial size T = 12 and rigid stationary 
ends bounding the annulus axially in order to compare with experiments For the main 
part (c.f. Sees. IIHI and IIV|) of this work we imposed, however, axially periodic boundary 
conditions at z = and z = V = 1.6. So the axial wavelength of the TVF and the SPI 
patterns investigated there is A = 1.6 and the wave number is k = 2n/X = 3.927. 

B. Numerical method 

The calculations were done with a finite differences method in the r — z plane combined 
with a spectral decomposition in (p. Since we have been studying also finite length cylinders, 
say, with lids bounding the annulus vertically, we do not use here a spectral decomposition in 
axial direction which for axially periodic systems is a generic alternative. The discretization 
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lm) has been done on staggered grids in 
ig| ] . It yields simple expressions for the 



(a FTCS - Forward Time, Centered Space algorit 
the r — z plane following the procedure of Ref. 
derivatives, it does not require boundary conditions for the pressure, and it avoids difficulties 
with boundary conditions for more than one velocity field component at the same position. 
We used homogeneous grids with discretization lengths Ar = Az = 0.05 which have showed 
to be more accurate than non-homogeneous grids. Time steps were At < 1/3600. 
Azimuthally all fields / — u,v,w,p were expanded as 

f(r,v,z,t)= Ur,z,t)e m *. (2.4) 

For the flows investigated here a truncation of the above Fourier expansion at m max = 8 
was sufficient to properly resolve the anharmonicities in the fields. The system of coupled 
equations for the amplitudes f m (r, z, t) of the azimuthal normal modes — m max < m < m max 
is solved with the FTCS algorithm. Pressure and velocity fields are iteratively adjusted to 
each other with the method of 'artificial compressibility' |20J] 

dpW = -/3V-uW (0</?<l) (2.5) 

p (n+l) = p (n) + dp (n) ^ 



= u (n) _ AtV{dp {n) ) . (2.7) 



The pressure correction dp^ in the n-th iteration step being proportional to the divergence 
of is used to adapt the velocity field vS n+1 \ The iteration loop (|2.5H2.7|) is executed for 
each azimuthal Fourier mode separately. It is iterated until V ■ u has become sufficiently 
small for each m mode considered - the magnitude of the total divergence never exceeded 
0.02 and typically it was much smaller. After that the next FTCS time step was executed. 

For code validation we compared SPI solutions with experiments and TVF solutions 
with previous numerical simulations |n| and close to onset also with Ginzburg-Landau 



results 



21(. Furthermore, we compared bifurcation thresholds of the nonlinear SPI and TVF 



solutions with the respective stability boundaries of the linearized NSE |6j, |14 1 obtained by 
a shooting method that is described in detail in Q| . As expected from our experience with 
primary vortex structures in the Taylor- Couette and Rayleigh-Benard problem lie the MAC 
FTCS bifurcation thresholds for our discretization typically 1 - 2 % below the respective 
linear stability thresholds. This deviation significantly reduces for finer discretizations. We 
also investigated how the nonlinear solutions change when varying m max and/or the grid 
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spacing. From these analyses we conservatively conclude that typical SPI frequencies have 
an error of less than about 0.2% and that typical velocity field amplitudes can be off by 
about 3 - 4%. Time steps were always well below the von Neumann stability criterion and 
by more than a factor of three below the Courant-Friederichs-Lewy criterion. 

In order to trace out the unstable parts of bifurcation branches of TVF and SPI solutions 
we applied different stabilization methods that are described in Sec. IIII B II 



III. SPIRAL VORTICES AND TAYLOR VORTICES 

In this section we first briefly review spatio-temporal properties of spiral vortices (M ^ 0) 
and Taylor vortices (M = 0) in the absence of any externally enforced axial through-flow. 
Here M is the azimuthal wave number of the respective vortex structure. Then we present 
our results on the bifurcation behavior of M = and M — ±1 vortex solutions and on their 
flow structure. 

They both grow out of the basic CCF state, vlccf = vccF(r)e v , that is rotationally 
symmetric, axially homogeneous, and time translationally invariant. Here in our system 
with 1] = 1/2 the radial profile of its azimuthal velocity reads 

VccF{r) = r + -. (3.1) 

A. Spatio-temporal structure 

The spiral vortex structure is periodic in tp, z, and t. It rotates uniformly as a whole 
like a rigid object in azimuthal direction thereby translating with constant phase velocity in 
axial direction — the spiral fields f(r,(p,z,t) do not depend on <p, z,t separately but only 
the phase combination 

<j> = kz + M<p-w(k,M)t. (3.2) 

Here k is the axial wave number that we always take to be positive and u(k, M) is the fre- 
quency. Thus, with f(r, if, z, t) = F(r, 0), the spiral pattern is one dimensional. Comparing 
the Fourier decompositions 

f(r, <p, z,t) = J2 fm,n(r, t) e^ + "^ = ^ F u (r) e** = F(r, 0) (3.3a) 

m,n v 
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one finds that 

f m ,n{r, t) = 5 mjnM e- in ^F n (r) . (3.3b) 

Thus only the mode combinations m = nM appear in a SPI with azimuthal wave number 
M. 

The SPI phase is constant, </>o, on a cylindrical surface, r = const, along lines given by 
the equation 

M u(k,M) 1 , > \ 

Thus, on the <p — z plane of such an 'unrolled' cylindrical surface these lines of constant 
phase are straight with slope —M/k. And an azimuthal wave number M > implies a left 
handed spiral (L-SPI) while M < refer to right handed spirals (R-SPI) with our convention 
of taking k to be positive. L-SPI and R-SPI being mirror images of each other under the 
operation z — > —z are symmetry degenerate flow states. Which of them is realized in a 
particular experimental or numerical setup depends on the initial conditions. 

The lines of constant phase and with it the whole spiral structure rotates in if with 
angular velocity 

<Pspi = ■ (3.5) 

Its direction strongly depends on the inner cylinder's rotation due to the influence of the 
CCF. The latter decisively determines the shape of the linear spiral eigenmodes that can 
grow beyond the stability boundary of the CCF state against perturbations with azimuthal 
wave number M^O. In the parameter range explored here the spirals rotate into the same 
direction as the inner cylinder, i.e., into positive ^-direction so that u(k, M)/M is always 
positive, i.e., u = sign(M)\u\. From this rigid rotation one immediately infers from Eq. ()3.2|) 
that the axial phase velocity 

w ph = — = -jt^Pspi (3.6) 

of an L-SPI (M > 0) is positive and of an R-SPI (M < 0) it is negative. 

For the rotationally symmetric (M = 0) structure of toroidally closed Taylor vortices 
the lines of constant phases are parallel to e v . This M = pattern is stationary (u — 0). 
The main reason is that the azimuthal flow of the basic CCF state being precisely parallel 
to the vortex lines of constant phase cannot advect them. However, an axial mean flow, 
being perpendicular to them can advect them: an externally enforced axial through-flow 
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of strength Re causes a non-zero axial phase velocity of the Taylor vortex pattern that 
grows linearly with Re, at least when phase pinning effects are absent as for axially periodic 
boundary conditions. 



B. Bifurcation behavior 



In the parameter regime considered here the bifurcation thresholds for nonlinear SPI and 
TVF solutions, i.e., the linear stability boundaries of the CCF state against M = ±1 and 
M = vortex perturbations |fj| differ only slightly from each other. For our fixed wave 
number of k = 3.927 they intersect at (Rf = 95.25, R\ = —73.69) where these two different 
vortex modes are bi-" critical" in the sense that their growth rates are simultaneously zero. 
The stability boundaries were obtained with a shooting method from the linearized NSE. The 
nonlinear SPI and TVF solutions that were determined with the numerical method described 
in Sec. Ill BI give bifurcation thresholds that differ as a result of the FTCS discretization 
errors by at most 2% from the linear stability analysis. However, this difference can grow 
with externally applied through-flow up to, say, 5% at Re ~ 40 (c.f. Sec. IIV|) when the 
discretization is not refined. 



1. Radial flow amplitudes of TVF and SPI 

The bifurcation of both, TVF and SPI solutions is forward as shown by the bifurcation 
surface over the R1 — R2 plane of Fig.|2j There the respective vortex solution is characterized 
by the primary Fourier amplitude, |w m , n |, of the radial flow intensity at mid gap, r = ri + 0.5, 
taken as order parameter with m denoting the azimuthal mode index and n referring to the 
axial one, respectively. Thus, Fig. |21 shows |ito,i| for the TVF solution by thin lines and 
\ u x,i\ = \ u -i,x\ f° r the two symmetry degenerate M = ±1 SPI solutions by thick lines, 
respectively. In each case stable (unstable) solutions are represented by full (dashed) lines. 
The different stability regions labelled A - E are explained in the caption of Fig. |2l 

The stability of the vortex states refers to our system with fixed axial periodicity length. 
Thus, e.g., Eckhaus or Benjamin- Feir instabilities [2^ that can destabilize periodic patterns 
in infinite and large systems do not occur here. Furthermore, our periodic boundary con- 
ditions allowing free phase propagation enhance the existence range as well as the stability 
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range of SPI solutions in comparison with, say, Ekman vortex generating stationary lids 
that axially close the annulus in an experimental setup. The latter suppress phase propa- 
gation in their vicinity so that phase generating and phase destroying defects near opposite 
boundaries are necessary for the realization of spirals in the bulk of such systems. 

In our setup TVF is for R 2 > R 2 stable close to onset. And it remains so at least up to 
the largest value of R\ = 130 shown in Fig. |21 for larger Ri TVF eventually undergoes an 
oscillatory instability. For more negative R2 < R 2 TVF is unstable at onset (region C in 
Fig. |2J) but becomes stable at larger Ri in region E. The unstable TVF solution branch in 
region C was obtained by suppressing any m^O modes in the field representation (|2.4jl . i.e., 
by allowing only rotationally symmetric solutions. Lifting this mode restriction infinitesimal 
m 7^ perturbations drive the system in the parameter region C of Fig. 121 away from the 
unstable TVF solution into a stable SPI state. 

Spirals, on the other hand, are for R 2 < R 2 stable close to onset and remain so at least 
up to the largest value of -Ri = 130 shown in Fig. |2] while for R 2 > R 2 they are unstable 
at onset (region D in Fig. |2J). But then they become stable at larger R\ in region E. The 
unstable SPI solution branch in region D was obtained by suppressing m = contributions 
to the radial velocity field u at mid gap location. This stabilized the SPI solution against 
the growth of TVF. Lifting this restriction of the available mode space the unstable SPI 
solutions in region D decay into stable Taylor vortices. 

In the relatively large region E both, SPI as well as TVF solutions coexist bistably and 
the final vortex structure to be found here depends on the initial conditions and the driving 
history of R±,R 2 . Note in particular that for our periodic boundary conditions the region E 
with stable spirals extends to positive R 2 , i.e., to a situation with co-rotating cylinders. 

2. SPI frequencies 

In Fig. El the spiral frequencies u are plotted over the same control parameter range as 
the radial flow amplitudes in Fig. |21 Also here we include - for the sake of comparison 
with Fig. |21 the identification of the different stability regions of TVF and SPI solution by 
the symbols A-E explained in the caption of Fig. |2l At onset u agrees within the numerical 
accuracy of our nonlinear code with the eigenvalue resulting from the linear stability analysis 
of the CCF state. 
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The nonlinear SPI frequencies further away from onset vary smoothly: the bifurcation 
surface of uj in Fig. |3] has the shape of a cloth that hangs down from a frame given by the 
linear onset spiral frequencies uj(Ri, s tab) at the stability threshold Rx, s tab{R2) of CCF. The 
location of minimal uj on the bifurcation surface is shown by a thick line in Fig. El Thus, 
the nonlinear SPI frequencies are typically smaller than the linear ones but do not deviate 
substantially from them. 

Since the linear onset frequencies show a characteristic variation along the bifurcation 
threshold, Rx^tabiRo) , that dictates the form of the whole uj bifurcation surface we discuss 
them in some detail. They, furthermore allow for a simple, yet semiquantitive explanation 
of the phenomenon of rigid body rotation of spirals in terms of a passive advection dynamics 
of M = ±1 vortex perturbations, e 1 ^, with lines of constant phase, = kz + Mip — cut, that 
are oriented obliquely to the "wind" of the basic azimuthal CCF. To that end we compare in 
Fig-IUthe onset spiral frequency uj{R\ t stab) a t the stability threshold R\, a tab{R2) of CCF with 
the "model" frequency uj mo< ^ e ^{Ri yS t a b) which is evaluated also at the stability threshold 
Ri,stab{R2)- Here 



is the mean of the rotation rate of the CCF, ujccf = vccf/t. For R 2 < the averaging is 
done over the radial domain between inner cylinder, ri, and the first zero, ro, of vccf{t) 
(jHHJ- Thus, at the stability threshold R\, s tab(R2) one has 



when R 2 < 0. However, when R 2 > 0, i.e., when vccf remains positive throughout the 
gap ro is replaced by r 2 - The restriction of the radial average to the range between r\ and 
r is motivated by an argument of largely hand-waving nature: the linear eigenfunctions 
for marginally stable SPI modes are somewhat centered to this range where the growth of 
vortex perturbations is supported. 

Fig. 0]shows that the onset spiral frequency oo(Ri >s tab) agrees perfectly well with the mean 
CCF rotation frequency (|H.7jl when R 2 > 0. For R 2 < the model ansatz ()3.7j) for the global 
spiral rotation rate overestimates slightly the spiral frequency since Eq. ()3.7j) does not contain 
contributions from negative CCF rotation rates between r and r 2 . In fact, if one extends 




(3.7) 



2R 2 — ARi :Stab 

2R 2 — R\ t stab 



(3.8) 
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in an ad hoc way the averaging domain slightly beyond r then the agreement improves 
significantly. Thus, the onset spiral frequency u(Ri :Sta b), i-e., the frequency eigenvalue can 
be seen as the mean rotation rate of the CCF - albeit weighted appropriately by the critical 
eigenfunctions. 

C. Flow structure of TVF and SPI 

In this section we elucidate the flow structure of spiral vortices in comparison with Taylor 
vortices. To that end we consider the radial velocity field. In Fig.EJwe show the axial profiles 
of u(z) at mid gap position for Ri = 130 being fixed and various R% that cover the whole 
interval between the bifurcation thresholds, c.f. Fig. Eland the inset of Fig. El Full (dashed) 
lines refer to negative (positive) i?2- in each case the axial position of maximal radial outflow 
is chosen to lie at z — 0.5A. For the sake of better visibility two axial periods of the vortex 
profiles are shown. 

1. Anharmonicity: TVF versus SPI 

Typically SPI are less anharmonic than TVF. Also the profiles of both are less anharmonic 
for positive R2 than for negative R2 and the degree of anharmonicity increases when R2 
becomes more negative. For the mirror symmetric TVF this anharmonicity growth comes 
from a widening (narrowing) of the axial range A in (A out ) of radial inflow over which u < 
(u > 0) and the corresponding decrease (increase) of the inflow (outflow) velocity. For the 
L-SPI that propagate in Fig. El into positive z-direction the anharmonicity grows mainly by 
flattening (steepening) of the wave profiles ahead of (behind) the crests. However, A in / A out 
increases also for SPI albeit less than for TVF. 

The variation of the anharmonicity of the vortex profiles can be read off more quantita- 
tively from the results of an axial Fourier analysis. To that end we show in Fig. El the ratios 
\u n /ui\ of the n-th and first axial Fourier modes of the profiles of Fig. El as a function of R2 
for fixed R\. With growing distances from the bifurcation thresholds at positive and nega- 
tive i?2 the anharmonicity grows for TVF as well as for SPI. It does so most precipitously 
near the thresholds at negative R2 of about -150 in Fig. El 

At negative R2 the anharmonicity of TVF can be for rapidly counter rotating cylinders 
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already close to threshold so large that \u2/ui\ > 1. This property reflects the fact that for 
sufficiently negative R2 Taylor vortices are effectively smaller in size than the gap width. 
There are two main reasons for this size reduction which are both connected to the tendency 
of vortices to have circular shape: (i) the axial periodicity length A = 1.6 reduces the axial 
vortex size relative to the gap and, more importantly, (m) the TVF intensity is radially 
restricted not to extend significantly beyond the zero of CCF at r since according to the 
Rayleigh criterion m = radial momentum transport is suppressed by opposite pressure 
gradients for r > r where the CCF stratification of the squared angular momentum density 
is stable. With R2 becoming more negative ro moves inwards and the radial size of Taylor 
vortices reduce. 

However, the m = Rayleigh criterion does not apply to SPI. Their m 7^ radial 
momentum transport extends further beyond r . Therefore SPI vortices fill out the whole 
gap more than Taylor vortices, c.f. Fig. and consequently they are less anharmonic. 

2. Mirror symmetry breaking of SPI 

TVF shows axial mirror symmetry around the position of maximal radial outflow, z = 
0.5A, in Fig. 03 In order to measure the degree to which this symmetry is broken in SPI we 
have used the asymmetry parameter 



evaluated at mid gap with z' = locating the largest radial SPI outflow at this r-value. In 
this way we found, e.g., for the spirals of Fig. that the smallest P ~ 0.2 occurs for spirals 
with the smallest frequency uo m i n ~ 23.4 at R2 — —74. Increasing R2 from this value all the 
way toward the upper SPI bifurcation threshold at R2 — 48 the frequency increases but P 
remains roughly unchanged at about 0.2. On the other hand, when decreasing R2 from -74 
the asymmetry parameter increases with increasing u up to P ~ 1 close to the lower SPI 
bifurcation threshold R2 — —158. Thus, fast propagating spirals at large negative R2 show 
the largest mirror symmetry breaking. 



P 



$\u{z')-u{-z')\dz' 
$\u(z')+u(-z r )\dz' 



(3.9) 
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D. Comparison with experimental results 

In order to check our numerical results we made a few comparisons with experiments. For 
exatnple, in F,g.®we 8 how the a^ia. proffle of the radial flow „ W of an L-SPI at + 0.4. 
Symbols denote Laser-Doppler velocimetry measurements |2| and the full line a numerical 
simulation, both done in a setup of height T = 12 with rigid, non-rotating lids at both ends 
of the annulus. In each case the spirals were monitored at mid-height of the cylinders where 
they had the common wavelength A ~ 1.76. Since absolute experimental velocities were not 
available we have scaled the experimental maximum in Fig. E to that of our simulation (full 
line). Without knowledge of the experimental error-bars we consider the agreement between 
symbols and full line to be satisfactory. 

The dashed line shows a numerical profile obtained for axially periodic boundary condi- 
tions imposing the wavelength A = 1.6. It differs slightly from the SPI profile (full line) in 
the bulk part of the T = 12 system with rigid ends. The difference is presumably related to 
the fact that the axial flow, and in particular the mean- flow wq (j4.4|l . is different in these 
two CHSGS clS discussed in Sec. II V A 21 

In Fig. El we compare the frequency variation of experimental and numerical L-SPI with 
R\. Symbols and the full line come from Laser-Doppler velocimetry measurements jj| and 
numerical simulations, respectively, of the aforementioned Taylor-Couette setup (77=0.5) of 
height T = 12 with rigid, non-rotating lids at both ends. Note that not only the frequency 
values of these experimental and numerical SPI states agree reasonably well with each other 
but also their existence range in R\. Its lower end marks the oscillatory onset. At the 
upper end in R\ these SPI lose their stability to TVF — in experiments as well as in the 
simulations. 

However, under axially periodic boundary conditions the existence range of stable SPI 
extends to significantly larger values of R\ lying outside of the plot range of Fig. El The 
dashed line in Fig. |H1 refers to simulations done with axially periodic conditions (A = 1.6) 
that allow for a free propagation of phase. And, in addition, they allow the Reynolds-stress- 
sustained mean axial flow wq ()4.4j) to have a finite net part < w > (j4.r>j) that is negative for 
our parameters - c.f. Sec. IIV A 21 In order to compare with the SPI frequencies for rigid end 
conditions we subtract from the oscillation frequencies under periodic boundary conditions 
(dashed line) the pure Galilean contribution (w)k and obtain the dash-dotted line. Note how 
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close the latter lies to the SPI frequencies in the system with rigid end conditions. Thus, we 
find that the SPI frequency differences Q] for the two different end boundary conditions 
are mostly due to whether the Galilean contribution (w)k is suppressed or not. 



IV. EXTERNAL THROUGH-FLOW 



Here we discuss the influence of an externally imposed axial through-flow on spiral and 
on Taylor vortices. Since the effect of an axial through-flow on TVF has been investigated 
for i?2 = in several works, we focus our investigation on SPI vortices. 

The through-flow is enforced by adding in the NSE for the axial velocity component a 
constant pressure gradient of size d z pAPF throughout the annulus. In the absence of any 
vortex flow, i.e., for sub-critical control parameters this pressure gradient, d z pAPF, drives an 
annular Poiseuille flow (APF) with a radial profile of the axial through-flow velocity given 
by 



/ v QzPAPF 
WAPFVr) = : 



' 2 , l+?7 1 , (1 + V) Ml - V) 1 
r + — In r + 



(4-1) 



[l — r])\nr] (1 — 77) In 7/ (1 — 7/) 2 

We checked that our numerical code reproduces this analytical solution (|4.1j) of the NSE. 
We use its mean to define the through-flow Reynolds number by 

/ / \ \ 7~> QzPapf 1 - rf + (1 + 7/ 2 )hi7y 

(w APF (r)} = Re = {1 _ r])Hn7] ■ (4-2) 

Hence positive (negative) Re implies an axial flow, wapf{t), in positive (negative) in- 
direction. The last equality in Eq. ()4.2|) establishes the relation between the externally 
applied additional axial pressure gradient and the through-flow Reynolds number Re. 



A. Counter-rotating cylinders 

Fig. 1101 shows how the through-flow influences L-SPI, R-SPI, and TVF at the fixed char- 
acteristic driving combination Ri = 120, R2 = —100 that is located in Figs. Eland El in the 
region C close to the border to region E. For this parameter combination TVF is unstable 
when Re = and it remains unstable in the .Re-range shown in Fig. El This is of relevance 
for the through- flow induced transitions between L-SPI and R-SPI (c.f. further below). 
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1. Bifurcation behavior 

We present in Fig- ITUT a) the primary Fourier amplitudes, |w m , n |, of the radial flow intensity 
at mid gap versus Re. These are |wi,i| for the M = 1 L-SPI, for the M = — 1 R-SPI, 

and |«o,i| f° r TVF. Fig. ITUTb) shows their axial phase velocity, w p h = co/k, and Fig. fTOlf c) 
shows the net mean axial flow 



For Re = the two spirals are mirror images of each other: their radial velocities are 
the same and all respective axial velocities have the same magnitude but opposite direction. 
Note that the SPI Reynolds stresses drive an axial flow to be discussed further below. Its 
net mean, (w) ()4.3j) . is directed opposite to the phase velocity, w p h, of the respective spiral 
when Re = 0. Note, however, the difference in size between \w p h\ — 7.1 and \(w)\ ~ 1.1 j^. 

A finite through-flow breaks the mirror symmetry between the M = 1 L-SPI and the 
M — — 1 R-SPI. Their radial flow amplitudes evolve with through- flow as shown in Fig. nffl a). 
We performed also a linear stability analysis of the combined CCF-APF state. It shows 
that for our control parameters R\ = 120, R2 = —100 the amplitudes of the M = ±1 SPI 
solutions go to zero at the bifurcation threshold values of Re = =f19.07 and Re = ±50.95. 
The numerical solutions of the full nonlinear NSE showed in addition that the L-SPI (R- 
SPI) is unstable near the first threshold, Re ~ —19 (Re ~ 19), and that it is stable near the 
second one, Re ~ 50 (Re ~ —50). 

For small through-flow - say, for —6 ^$ Re ^ 6 in Fig. - the two spiral solutions 
coexist bistably; their particular realization depending on initial conditions. However, with 
increasing |i?e| that spiral suffers a thro ugh- flow enforced loss of stability for which the 
phase velocity changes sign. This happens roughly when the through-flow has become 
sufficiently strong to revert an originally adverse axial phase propagation. For example, the 
M = —1 R-SPI of Fig. fTOl propagate for small Re < 6.6 axially downwards (i.e. opposite 
to the externally imposed through- flow) as for Re = 0, then become stationary, and finally 
propagate upwards in through-flow direction for Re ^ 6.6. Similarly, by symmetry, the 
M — 1 L-SPI propagates in a small negative through-flow upwards against the through-flow 
for Re ^ —6.6 and downwards, i.e., in through-flow direction for Re < —6.6. 

The direction of the imposed through-flow is the preferred one for stable phase prop- 




(4.3) 



n 
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agation: A spiral that has started at small |i?e| to move against the wind dies out - or, 
more precisely, becomes unstable - when the wind becomes sufficiently strong to turn it 
back. Only that SPI is stable at large \Re\ ^ 7.2 in Fig. Elthat keeps propagating into the 
preferred direction of the through-flow. The other one is unstable at large \Re\. 

The through-flow enforced loss of stability of one SPI state and the transition to the 
remaining stable one is indicated schematically in Fig. ITUT a) by vertical arrows. However, 
we should like to stress that the transition is somewhat complex extending over the through- 
flow interval 6 < \Re\ ^ 7.2 the center of which locates the zero of w p h at |i?e| ~ 6.6. In this 
interval there are stable, mixed states with finite L- and R-SPI modes. Their amplitudes 
seem to vary largely continuously with Re (with possibly some saddle-node discontinuity) 
between the pure SPI solutions: the amplitude of the spiral that loses the stability com- 
petition decreases with growing \Re\ towards zero while the amplitude of the winning one 
increases from zero to the pure monostable final SPI state. 

Note that since TVF is unstable for the parameters of Fig.Elit does not offer an alterna- 
tive transition to a final M = state as for the parameters of Sec. IIVBI There, for R2 = 0, 
the through-flow induces a transition to stable TVF rather than to the stably coexisting 
SPI with preferred propagation direction. Only when TVF is eliminated there does the 
transition occur to the then monostable spiral — for details see Sec. IIVBI 

We also made a few calculations in a regime where TVF stably coexists with SPI for 
counter-rotating cylinders. Also then the through-flow induces preferably a transition to 
stable TVF rather than to the stable SPI state. Thus, when the through-flow destabilizes, 
e. g., the M — — 1 R-SPI then typically the M = TVF modes grow rather than the M — 1 
L-SPI modes. 

2. Axial velocities w p h,WQ, and (w) 

In the through- flow range shown in Fig. ^] the phase velocity w p h and the net mean flow 
(w) vary roughly linearly with Re. The slopes dw p h/dRe and d{w)/dRe for SPI as well as 
for TVF are roughly 1. 

While the phase of the M = ±1 SPI reverts its propagation direction at Re ~ =f6.6 the 
net mean flow changes sign already at Re ~ ±1.2. The reversal of the latter does not seem 
to have any consequence. But the through-flow enforced reversal of the phase velocity seems 
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to be responsible for the destabilization of the SPI that propagate at small \Re\ against the 
wind, i.e., into the "wrong" direction. 

In Fig. ^2 we show how the radial profiles of the mean axial flow 



2tt 

1 



w (r) = — J w(r,(p,z,t)dcp (4.4) 
o 

of spirals shown in Fig. evolve with the through-flow in the range —4 < Re < 14. We 
checked that Wq is independent of z and t and that our spirals propagating in the externally 
imposed axial pressure gradient still have the SPI symmetry, i.e., the flow fields depend on 
z, <p, t only via the phase combination <fi ()3.2|) with an oscillation frequency ui that is modified 
by the through-flow. Then one finds from the NSE for the m = azimuthal mode of the 
axial velocity field, 

d r + -) d r w = ( d r + - J (uw) + d z p , (4.5) 



that the SPI mean flow can be driven by Reynolds stresses and/or by mean axial pressure 
gradients. For Re = the pressure is enforced to be axially periodic, hence d z po(Re = 0) = 0. 
So in that case the mean axial flow is driven solely by the nonlinear Reynolds stresses. They 
are rather large. For example for the R-SPI propagating at Re = in negative ^-direction 
with phase velocity w p h — —7.1 the maximum of w (r) is about 3, i.e, directed opposite to 
the phase propagation and almost half as large in magnitude as w p h- The net mean flow (w) 
()4.3|) is for this case still about 1.1 and also opposite to w p h- 

As an aside we mention that rigid axial end conditions enforce (w) = throughout the 
annulus. They generate an adverse axial pressure gradient that compensates the Reynolds 
stresses so that Wq is practically zero in the bulk part where SPI are realized. Only in 
the Ekman region wq becomes finite showing TVF behavior there. 

For the R-SPI of Fig. propagating at Re > opposite to the external through-flow 
the maximal mean flow is located roughly at mid-gap. However, for the SPI propagating 
into the direction of the external through-flow, i.e., the R-SPI for Re < and the L-SPI for 
Re > the extremum of wq(t) is shifted towards the inner cylinder. The mean flow profiles 
of the spirals of Fig. ^2 are given within about 5% by the superposition 

w (r; Re) ~ w (r; Re = 0) + w APF (r; Re) (4.6) 

of the pure, Reynolds stress generated flow w (Re = 0) of the respective SPI plus the pure, 
pressure gradient enforced APF flow wapf{.R^) (j4.1|) . This holds for L-SPI as well as for R- 
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SPI, irrespective of whether they propagate into the direction of the through-flow or against 
it. 

3. Spiral profiles 

The through-flow changes the structure of the SPI. This is documented in Figs. ^] and 
fTHl The arrows in Fig. II 21 representing the u, w vector field of L-SPI in the r — z plane show 
the effect of imposing an axial through-flow that increases from Re = —5(a) to Re = 10(d) 
in steps of 5. Note, however, that the externally imposed axial pressure gradient does not 
just add wapf(t) to the axial velocity field w. It also modifies all vector field components of 
the SPI. The axial profile of the radial flow u(z) for example is changed by the through-flow 
as shown in Fig. ^] for increasing Re. Here the axial asymmetry of the upwards propagating 
L-SPI is reduced by steepening up the leading part of u(z) ahead of the wave crests. This 
reduction of the mirror-asymmetry of the radial flow of L-SPI grows somewhat linearly with 
increasing Re. As an aside we mention that on the other hand the TVF profiles of u(z) 
become with increasing Re more and more asymmetric — the mirror asymmetry parameter 
P ()3.9|) increases for TVF linearly with Re. 

B. Non-rotating outer cylinder 

We have investigated the influence of an externally imposed axial through-flow on TVF 
and SPI also for stationary outer cylinder, R2 = 0. 

1. Bifurcation behavior 

In Fig.Elwe show the bifurcation behavior of TVF and SPI as a function of through-flow 
Reynolds number Re for R2 = 0, Ri = 100. This parameter combination lies well within 
the region E of Fig. El in which TVF, L-SPI, and R-SPI are all stable at Re = 0. 

Switching on the through-flow one sees in Fig. HUf a) how the dominant modes of these 
vortex structures vary with Re. That SPI loses its stability for which the through-flow 
enforces a reversal of the phase propagation as in the case of counter rotating cylinders 
(Fig. Unj). Thus, also here the direction of the imposed through-flow is the preferred one 
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for stable SPI at large \Re\. A spiral that has started at small |i?e| to move against the 
through-flow becomes unstable when the latter becomes sufficiently strong to turn it back. 
On the other hand that SPI remains stable at large |i?e| that keeps propagating into the 
preferred direction of the through-flow. 

As in Fig. El the loss of stability takes place in the vicinity of the Reynolds number 
where the axial phase velocity w p h [Fig. HHf b)] of the respective SPI goes through zero. This 
happens in Fig.E)f° r the M = ±1 SPI at Re ~ =f6.4. However here we found the transition 
from the then unstable SPI to occur to the stable TVF solution [c.f. arrows in Fig. Hffi a)] 
rather than to the other stable SPI. 

We have also investigated briefly the situation where the TVF solution was numerically 
eliminated (here, suppressing m = modes of the w-field at mid-gap position turned out 
to be an efficient way to globally reduce TVF towards zero). Also then the SPI that is 
unfavored by the through-flow loses its stability. However, with TVF being unavailable as 
final state the transition occurs in this case to the favored SPI in a way that seems to be 
similar to the one described in Sec. IIV A II 

Without the above described numerically imposed mode restriction TVF is stable for 



26, 
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For 



moderate through-flow rates while at sufficiently large \Re\ SPI are stable 
our parameters TVF decays at Re ~ ±34 into a M = ±1 SPI as indicated by arrows in 
Fig. da). 

For small through-flow the phase velocity w p h and the net mean flow (w) vary roughly 
linearly with Re. The initial slopes dw p h/dRe and d(w)/dRe are for SPI as well as for 
TVF roughly 1. However, at larger Re one sees in Fig. IT4T c) that in particular (w) shows 
nonlinear corrections. 



2. Phase diagram 

Fig.ll5lshows the phase diagram of TVF, R-SPI, and L-SPI for stationary outer cylinder in 
the control parameter plane spanned by Re and R\. The existence range of the vortex states 
is bounded from below by the bifurcation threshold (full line in Fig. ITSJl of the respective 
vortex solution out of the combined CCF-APF basic state. These bifurcation thresholds 
result from a linear stability analysis of the CCF-APF state The one for TVF increases 
quadratically for small Re. Also the SPI threshold curves in Fig. ED have a somewhat 
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parabolic shape, however with minima shifted to finite Re. Thus, the threshold for L-SPI 
first decreases for small positive Re but eventually increases at larger Re. By symmetry the 
R-SPI threshold curve in Fig.E3is the mirror image under Re — > —Re of the L-SPI threshold 
curve. Hence small through-flow destabilizes (stabilizes) the CCF-APF state against spirals 
that propagate into (against) the through-flow direction. 

Note that for small Re in Fig. TVF bifurcates first when increasing Ri. But for 
sufficiently large Re the bifurcation sequence of TVF and SPI is reversed since the bifurcation 
threshold for TVF curves up faster with increasing Re than the one for L-SPI. After their 
intersection stable SPI bifurcate first out of the CCF-APF state. Hence, for example in 
region E of Fig. ED only stable L-SPI exist, in region D TVF exists but only as unstable 
solution, and in region B they exist bistably. 

The dashed lines in Fig. are stability boundaries of the vortex solutions. Different re- 
gions of Fig. bet ween various stability boundaries and bifurcation thresholds are identified 
with the respective stability properties of the vortex states in the caption of Fig. [T^l 

V. SUMMARY 

We have numerically simulated vortex flow structures of different azimuthal wave numbers 
M in the Taylor-Couette system with counter-rotating as well as with co-rotating cylinders. 
In particular we have investigated the effect of an externally imposed axial through-flow on 
the spatio-temporal properties and on the bifurcation behavior of M = 1 L-spirals, M = — 1 
R-spirals, and M = Taylor vortices. 

To that end we first have determined for zero through-flow, Re = 0, the bifurcation sur- 
faces of the appropriate order parameters characterizing SPI and TVF solutions over the 
Ri — R2 control parameter plane of the inner and outer cylinder's Reynolds numbers. For 
the parameter combinations explored in this work these bifurcations out of the basic CCF 
state are forward and their order of appearance determines the stability of the respective 
bifurcating vortex state: the vortex solution that bifurcates second is unstable. But it even- 
tually becomes stable with increasing distance from the bifurcation threshold so that, e.g., 
for larger Ri there is a large region in the R\ — R2 plane with bistability of TVF and SPI. In 
particular the existence region of stable SPI extends for axially periodic boundary conditions 
even to positive R2 with co-rotating cylinders. Unstable solution branches were obtained by 
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selectively suppressing destabilizing modes. Stable ribbons, i.e., nonlinear combinations of 
M — ±1 spirals were not found. 

Simulations of axially finite systems with rigid, non-rotating lids showed in good agree- 
ment with experiments how the stable existence range of SPI is reduced by stationary Ekman 
vortices which suppress phase propagation at the two ends. Also the frequencies and the 
wave profiles of the spiral vortices in the bulk of the numerical and experimental systems 
agreed well with each other. Spiral profiles obtained for periodic and rigid end conditions 
do not differ much. On the other hand, the respective frequencies differ basically by the 
Galilean contribution (w)k. Here (w) is the net axial mean flow that the nonlinear Reynolds 
stresses of a spiral with axial wave number k sustains with axially periodic end conditions 
but not with impermeable ends. 

Furthermore, we showed how the phenomenon of rigid body rotation of spirals can be 
understood quantitatively in terms of the passive advection dynamics of M — ±1 vortex 
perturbations whose lines of constant phase are oriented obliquely to the azimuthal CCF. The 
onset spiral frequency is the mean rotation rate of the CCF, albeit weighted appropriately 
by the critical eigenfunctions with the consequence that L-SPI as well as R-SPI rotate into 
the same direction as the inner cylinder. The nonlinear SPI frequencies are typically smaller 
than the linear ones but do not deviate substantially from them. 

A finite through-flow breaks the mirror symmetry between the L-SPI and the R-SPI and 
it changes the structure of the SPI. The externally imposed axial pressure gradient does not 
just add the annular Poiseuille flow wapf(^) to the axial velocity field. It modifies the SPI 
structure, e.g., the profiles of the radial flow in a characteristic way. 

For Re = L-SPI propagate axially upwards and R-SPI downwards. When they are 
initially stable they continue to coexist bistably for small through-flow. However, they are 
no longer mirror images of each other and their phase velocities differ by an amount oc Re. 
Then, with increasing \Re\ that spiral loses its stability for which the through- flow enforces 
the phase velocity to change direction. Only that SPI is stable at large \Re\ that keeps 
propagating into the preferred direction of the through-flow. The other one is unstable at 
large \Re\. 

The SPI that loses stability upon reverting its propagation direction — i.e. the R-SPI 
(L-SPI) for positive (negative) Re — preferentially undergoes a transition to propagating 
TVF provided the latter is available as stable vortex state. Otherwise the transition is to 
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the then monostable L-SPI (R-SPI). Such a situation was explored in detail for negative 
R 2 where TVF was unstable and for other parameter combinations where the TVF solution 
was eliminated numerically. 

Also the situation where initially at Re = all three vortex solutions are stable was elu- 
cidated for different R\ — R 2 parameter combinations and in more detail for stationary outer 
cylinder, R 2 = 0. Here, a complete phase diagram was determined in the control parameter 
plane spanned by Re and R x . We found that small through-flow destabilizes (stabilizes) the 
basic CCF-APF state against spirals that propagate into (against) the through-flow direc- 
tion. For sufficiently large Re the bifurcation sequence of TVF and SPI is reversed since the 
bifurcation threshold for TVF curves up faster with increasing Re than the one for L-SPI. 
After their intersection stable SPI bifurcate first out of the CCF-APF state. Then there 
opens up a region at sufficiently large positive Re in which only stable L-SPI but no Taylor 
vortices exist for stationary outer cylinder. 
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R-SPI (M=-l) 




FIG. 1: Lines of constant phases, <j> = const, for spirals in the <p — z plane. Arrows indicate their 
velocities. 
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FIG. 2: Order parameter bifurcation surfaces of TVF (thin lines) and of SPI (thick lines) over 
the Ri — i?2 plane. Shown are primary Fourier amplitudes, |itf7in|> of the radial flow intensity at 
mid gap, r = n + 0.5, with axial mode index n = ±1. The azimuthal one is m = for TVF and 
m = ±1 for SPI, respectively. In each case full (dashed) lines denote stable (unstable) solutions. 
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FIG. 3: Bifurcation diagram of M = ±1 spiral frequencies oj over the R\ — R2 plane. The thick line 
locates the minima. The different stability regions A -E of TVF and of SPI solutions (c.f. caption 
of Fig. |2J) in the i?i — R2 plane are included for better comparison with Fig. |21 




FIG. 4: Linear frequency oj(Ri )S tab) of M = 1 spiral at onset, Ri^tabi^), in comparison with 
frequency u> model (Ri yS t a b) Q3.7[) resulting from rigid-body rotation model. 
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FIG. 5: Axial profiles of the radial velocity u(z) at mid gap position for R\ = 130 and various R2 
(along the thick horizontal line in the inset) covering the whole interval between the bifurcation 
thresholds marked TVF and SPI, respectively, in the inset; see also Fig. El Full (dashed) lines refer 
to negative (positive) i?2- In each case the maximal radial outflow is chosen to lie at z = 0.5A. 
For better visibility two axial periods of the vortex profiles are shown. The M = 1 L-SPI are 
propagating in positive z-direction. Parameters are 77 = 0.5, k = 3.927. 
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FIG. 6: Anharmonicity of TVF and SPI. The ratios |u n /^i| of the axial Fourier modes of the 
profiles of u(z) shown in Fig. are displayed here as functions of R2 for fixed R\ = 130. The 
bifurcation thresholds are located at the zeroes. Parameters are rj = 0.5, k = 3.927. 
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FIG. 7: Velocity field («, w) of TVF (left) and L-SPI (right) in an r - z plane. Vertical lines 
locate the zero of the azimuthal CCF flow vccf(i')- Parameters are r\ = 0.5, k = 3.927, R\ = 120, 
R 2 = -100. 
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FIG. 8: Comparison of experimental and numerical axial profiles of the radial velocity u{r\ + 0.4, z) 
of an L-SPI. For better visibility more than one period is shown. Symbols and the full line denote 
Laser-Doppler velocimetry measurements 9] and numerical simulations, respectively, of a Taylor- 
Couette setup of height T = 12 with rigid, non-rotating lids at both ends. Both refer to the 
bulk region at mid height with a common local wavelength of A ~ 1.76. There the experimental 
maximum of u is scaled to our simulation result. Dashed line refers to a simulation done with 
axially periodic conditions imposing a wavelength of A = 1.6. Common parameters are r\ = 0.5, 
R± = 111 with i?2 = —95 for the experiments and R2 = —96 for the simulations. 
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FIG. 9: Comparison of the frequency variation of experimental and numerical L-SPI with R\. 
Symbols and the full line come from Laser-Doppler velocimetry measurements ij] and numerical 
simulations, respectively, of a Taylor-Couette setup (77=0.5) of height T = 12 with rigid, non- 
rotating lids at both ends that enforce the net mean axial flow (w) (|4.3|) to vanish. Dashed line 
refers to a simulation done with axially periodic conditions (A = 1.6). They allow for a finite 
Reynolds-stress-sustained {w) that is negative for our parameters. Upon subtracting this Galilean 
contribution (w)k from the oscillation frequency under periodic boundary conditions (dashed line) 
one obtains the dash-dotted line that lies close to the SPI frequencies with rigid end conditions. 
Common parameters are R2 = —96, however, R2 = —100 for the full line. 




FIG. 10: Influence of an external through-flow on vortex structures, (a) Primary Fourier amplitudes 
of the radial flow field at mid gap for the M = 1 L-SPI (u± i), the M = —1 R-SPI i), and for 
TVF («o,i)- (b) Axial phase velocity uip/i = u/k. (c) Net mean axial flow (u>) (14.3(1 . Full (dashed) 
lines with filled (open) symbols refer to stable (unstable) states. Arrows indicate transitions after 
loss of stability, see text for details. TVF is unstable in the -Re-range shown here for our parameters 
Ri = 120, R 2 = -100, rj = 0.5, k = 3.927. 
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FIG. 11: Radial profiles of the axial mean flow ^o( r ) Q4.4JI of spirals shown in Fig. ^Jfor axial 
Reynolds numbers — 4 < Re < 14 increasing in steps of 2. Thick line refers to Re = 0. The 
transition from R- to L-SPI occurs around Re ~ 7, c.f. text. Parameters are -Ri = 120, R2 = 
-100,7] = 0.5, k = 3.927. 
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FIG. 12: Velocity field (u, w) of L-SPI in an r — z plane for Re=-5 (a), 0(b), 5(c), 10(d) . Parameters 
are rj = 0.5, k = 3.927, i?i = 120, R 2 = -100. 
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FIG. 13: The effect of an external through- flow on the axial profiles of the radial velocity of L-SPI. 
Lines show u(z) at mid gap position for Re = —5 to Re = 20 in steps of 5. Thick one refers to 
Re = 0. In each case the maximal radial outflow is chosen to lie at z = 0.5A. Parameters are 
Ri = 130, R 2 = -100, 7] = 0.5, A = 1.6. 
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FIG. 14: Influence of an external through-flow on vortex structures, (a) Primary Fourier amplitudes 
of the radial flow field at mid gap for the M = 1 L-SPI (mi 1), the M = — 1 R-SPI (u-i t \), and for 
TVF (uo,i)- (b) Axial phase velocity = u/k. (c) Net mean axial flow (w;) — i?e. Full (dashed) 
lines with filled (open) symbols refer to stable (unstable) states. Arrows indicate transitions after 
loss of stability, see text for details. Parameters are R\ = 100, R2 = 0, rj = 0.5, k = 3.927. 




FIG. 15: R\ — Re phase diagram of TVF, R-SPI, and L-SPI for stationary outer cylinder. Solid lines 
represent linear stability thresholds of the basic state, i.e., bifurcation thresholds of the respective 
vortex solutions out of the combined CCF-APF. Dashed lines are stability boundaries of the vortex 
states. The phase diagram is symmetric under Re — > —Re. Parameters are R2 = 0, rj = 0.5, k = 
3.927. 
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s: stable; u: unstable; -: nonexistent. 



